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ABSTRACT 


In  order  to  investigate  the  effects  of  gear  train  backlash  on  the 
performance  of  second  order  servomechanisms,  an  automatic  digital  com¬ 
puter  was  utilized  to  obtain  solutions  of  numerous  second  order  servo 
problems.  Solutions  were  obtained  for  variations  of  motor  and  load 
inertias,  motor  and  load  frictions,  combined  system  damping  coefficient, 
and  backlash  angle,  throughout  the  practical  range  of  these  parameters. 
The  data  generated  by  the  digital  computer  solution  has  been  analyzed 
and  presented  in  graphical  form  for  u?e  in  estimating  the  performance 
of  an  arbitrary  second  order  servo  system  with  backlash. 

The  program  via s  conducted  on  the  Control  Data  Corporation  Model 
1604  Computer  located  at  the  U.  S.  Naval  Postgraduate  School,  Monterey, 
California. 

The  writers  wish  to  express  their  appreciation  for  the  assistance 
and  encouragement  given  them  by  Professor  George  J.  Thaler  of  the  U.  S. 
Naval  Postgraduate  School. 
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Chapter  I 


INTRODUCTION 

A  servomechanism  which  is  designed  using  linear  theory  rarely 
performs  as  estimated.  The  discrepancies  between  design  performance 
and  actual  performance  are  caused,  in  general,  by  certain  non-line¬ 
arities  which  are  inherent  in  all  systems.  Of  the  many  non-line¬ 
arities  which  are  found  in  servo  systems,  the  one  under  investigation 
here  is  backlash. 

In  a  servo  system,  whenever  torques  are  transmitted  from  a  motor 
to  a  load  through  a  gear  train,  there  will  exist  a  certain  amount  of 
backlash  between  the  gears  as  shown  on  figure  1.  Present  day  gear 
drives  cannot  transmit  loads  continuously  without  some  backlash.  If 
the  gears  are  too  closely  meshed  they  might  land,  wear  excessively,  or 
introduce  random  torque  variations  appearing  as  noise.  This  is  partic¬ 
ularly  true  of  gear  trains  which  operate  over  wide  temperature  varia¬ 
tions.  Backlash  is  often  controlled  so  as  to  introduce  a  slight  dither 
to  a  system  to  provide  good  lubrication  to  overcome  stictional  effects. 

It  is  generally  known  that  backlash  in  a  servo  system  is 
destabilizing.  However,  the  estimation  of  the  effect  of  backlash  on 
the  transient  response  and  the  ultimate  steady  state  operation  of  a 
system  presents  a  very  difficult  problem.  As  soon  as  backlash  is 
introduced  into  an  otherwise  linear  system,  response  is  then  dependent 
upon  the  division  of  inertia  between  motor  and  load,  the  division  of 
viscous  friction  between  motor  and  load,  the  damping  coefficient  of  the 
system  when  operating  as  a  unit,  the  amount  of  energy  dissipated  as 
heat  when  contact  is  made  between  the  gear  teeth,  and  other  phenomena 
which  cannot  be  predicted.  Since  the  system  performance  is  affected 
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differently  in  the  presence  of  backlash  when  a  different  set  of  variables 
is  used,  each  system  presents  a  new  problem. 

In  practice,  knowledge  of  the  performance  of  a  system  may  be  desir¬ 
able  prior  to  the  construction  of  a  mock-up  or  simulation  system.  In 
order  to  obviate  the  necessity  of  estimating  the  performance  by  labori¬ 
ous  and  complicated  analytical  techniques,  design  charts  are  presented 
in  this  thesis  from  which  the  design  engineer  can  predict  the  magnitude 
of  maximum  overshoot  for  any  arbitrary  second  order  servo  with  backlash. 

A  qualitative  analysis  of  the  effect  of  backlash  on  the  time  of  maximum 
overshoot  and  the  settling  time  is  also  presented^ 

Previous  work  has  been  done  in  this  area  at  the  Naval  Postgraduate 
School  by  W.  J.  Lutkenhouse,  (Ref.  1),  and  N.  C.  New,  (Ref.  2). 
Lutkenhouse  solved  the  problem  by  construction  of  the  phase  trajectory 
on  the  phase  plane.  New  extended  the  work  of  Lutkenhouse  by  employing 
an  automatic  digital  computer  for  the  solution  of  the  problem,  and  was 
able  to  collect  a  large  volume  of  data  for  use  in  the  construction  of 
design  charts  for  predicting  the  size  of  limit  cycles  in  the  system. 

The  work  done  by  Lutkenhouse  and  New  was  limited  to  the  existence 
and  size  of  limit  cycles  in  a  system.  This  thesis  extends  that  work 
to  include  transient  behavior  of  load  motion.  The  computer  program 
written  by  New  was  utilized  for  the  solution  of  a  large  number  of 
problems,  the  data  from  which  was  subsequently  analyzed  and  processed 
for  the  construction  of  design  curves. 
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Chapter  II 
BACKGROUND 

Prior  to  the  development  of  high  speed  digital  computers,  it  was 
necessary  for  the  engineer  to  use  long  and  laborious  methods  for 
analyzing  the  behavior  of  servomechanism  systems  containing  non-linear 
elements.  Previous  investigation  of  methods  used  to  analyze  non-linear 
systems  has  resulted  in  several  approaches. 

One  of  the  most  widely  used  methods  for  analyzing  the  behavior  of 
a  non-linearity  is  the  phase  plane  diagram.  At  present,  this  method 
is  limited  to  control  systems  represented  by  a  second  order  differential 
equation  in  which  the  coefficients  can  be  constants  or  functions  of  the 
dependent  variable.  The  solution  to  this  problem  is  one  in  which  time 
is  not  expressed  explicitly,  but  only  the  velocity  and  position  of  the 
dependent  variable  for  each  instant  of  time  are  showp.  From  Ref.  1, 
it  can  be  noted  how  tedious,  time  consuming,  and  often  inaccurate  such 
a  method  is  for  obtaining  the  transient  response  of  a  second  order 
servo  with  some  non-linearity. 

Another  method  widely  used  in  non-linear  problems  is  the  use  of 
describing  functions.  This  method  uses  the  concept  of  obtaining  an 
"equivalent"  transfer  function.  If  a  sinusoidal  input  of  constant 
amplitude  is  applied  to  a  non-linear  system,  the  output  will  contain 
a  number  of  harmonics  of  the  input  frequency.  If  the  fundamental 
component  of  the  output  is  assumed  to  be  of  the  same  frequency  as  the 
input,  then  an  "equivalent"  transfer  function  can  be  obtained  which 
is  the  ratio  of  the  fundamental  component  of  the  output  to  the  sinusoi¬ 
dal  wave  input.  Once  this  relationship  of  the  output  to  the  input 
of  a  non-linear  system  is  obtained,  it  is  then  used  in  the  same 
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manner  as  the  transfer  function  in  the  linear  system.  Since  this 
method  assumes  that  the  higher  harmonics  have  negligible  effect  on  the 
system  and  that  the  frequency  of  the  input  is  the  same  as  the  output, 
it  is  seen  that  it  can  only  be  used  for  a  rough  approximation  in 
determining  the  amplitude  and  period  of  the  output.  This  method  is 
also  disadvantageous  because  of  the  difficulty  in  the  computations 
associated  with  the  determination  of  the  describing  function  for  a 
wide  range  of  amplitudes  and  frequencies.  There  is  no  simple  method 
for  evaluating  the  accuracy  of  a  describing  function  and  no  good  way 
to  check  the  validity  of  solutions  obtained  with  it. 

With  the  advent  of  the  analogue  computer,  another  method  of 
calculating  system  transient  performance  was  available.  This  method 
merely  simulates  the  problem  of  the  servo  and  its  non-linearities  to 
obtain  the  response  of  the  system.  This  method  is  of  good  value  in 
the  preliminary  evaluation  of  the  system  behavior  where  accuracy  is 
not  too  important.  However,  due  to  the  inherent  inaccuracies  of  the 
analogue  computer  this  method  does  not  give  precise  results. 

By  utilizing  the  digital  computer  in  the  solution  of  non-linear 
servos^ it  is  possible  to  overcome  many  of  the  disadvantages  inherent 
in  the  other  methods.  A  non-linear  system  which  can  be  described  by 
a  few  accurate  or  approximate  differential  equations  can  be  solved 
on  a  digital  computer  resulting  in  transient  and  steady  state  informa¬ 
tion.  In  this  thesis  a  digital  computer  was  used  to  solve  many  second 
order  servo  problems  with  backlash  present.  From  the  wealth  of  results 
obtained,  the  response  of  any  arbitrary  second  order  servo  with  back¬ 
lash  may  be  predicted.  It  is  believed  that  this  procedure  could  be 
applied  to  servos  of  higher  order  and  with  other  non-linearities. 
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Chapter  III 

DISCUSSION  OF  SYSTEM  BEHAVIOR 


The  system  'under  investigation  in  this  thesis  is  a  simple  second 
order  position  servo  control  having  the  motor  coupled  to  the  load  through 
gear  trains  having  backlash.  The  configuration  of  the  system  is  illus¬ 
trated  in  block  diagram  form  in  figure  2.  Position  feedback  is  taken 
from  the  load  shaft  including  the  backlash  within  the  feedback  loop. 

When  backlash  is  outside  the  feedback  loop,  i.e.,  position  feedback 
from  the  motor  shaft,  the  system  is  stable  and  does  not  result  in  a 
limit  cycle;  however,  it  is  subject  to  a  residual  steady  state  error 
in  response  to  a  step  input  which  error  may  be  as  great  as  the  magnitude 
of  the  backlash.  When  included  in  the  feedback  loop,  backlash  is  de¬ 
stabilizing  and  the  effect  is  dependent  upon  many  variables.  When  a 
step  position  input  is  imposed  on  the  system  there  is  initially  a  period 
when  the  backlash  is  taken  up,  and  the  system  performs  as  a  linear  sys¬ 
tem.  If  the  linear  system  damping  coefficient  is  less  than  the  critical 
value,  the  error  signal  will  eventually  cause  the  motor  shaft  to  reverse 
its  direction  of  rotation.  A  separation  between  the  motor  and  the  load 
will  occur  in  the  cycle  of  motion  whenever  the  deceleration  of  the  motor 
becomes  greater  than  that  of  the  load.  However,  if  the  frictional 
forces  on  the  load  produce  a  deceleration  greater  than  the  maximum  de¬ 
celeration  of  the  motor,  the  motor  and  the  load  will  maintain  '•ontact 
during  this  period  and  become  separated  only  when  the  motor  velocity 
reverses.  The  precise  instant  at  which  motor  and  load  separate  depends 
upon  the  nature  of  the  frictional  forces  acting  and  the  division  of 
inertia  between  the  motor  and  load.  When  the  gears  do  break  contact, 
the  load  drifts  and  the  system  is  no  longer  a  closed  loop  system  but  is 
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driven  open  loop.  Whenever  the  total  amount  of  backlash  is  taken  up 
and  gear  contact  is  made,  the  action  which  follows  is  very  complex. 

The  basic  laws  which  govern  the  ensuing  motion  are  the  law  of  conserva¬ 
tion  of  energy  and  the  law  of  conservation  of  momentum.  In  this  thesis 
the  assumption  is  made  that  the  contact  is  a  plastic  one  ’with  no  gear 
bounce,  and  that  energy  is  dissipated  as  heat  due  to  the  plastic  contact. 
Thus  only  the  law  of  conservation  of  momentum  are  applied.'*'  Following 
the  gear  contact  phase,  the  system  is  again  combined,  and  the  cycle  of 
events  begins  again. 

The  motion  of  the  motor  and  load  is  illustrated  in  figures  3  and  4. 
Angular  velocity  versus  angular  position  is  plotted  and  the  phase  tra¬ 
jectories  for  both  motor  and  load  are  shown  for  two  conditions  of  inertia 
and  friction  ratios.  In  figure  3  -it  is  noted  that  at  point  A,  the  sys¬ 
tem  reacts  to  a  step  input  of  1  radian  by  moving  in  a  trajectory  in 
accordance  with  the  differential  equation  for  the  linear  system.  At 

point  B  the  gears  which  couple  the  motor  to  the  load  separate;  the 

dotted  trajectory  from  B  to  C'  showing  the  path  followed  by  the  motor 
and  the  solid  line  from  B  to  C  showing  the  path  of  the  drifting  load. 

The  gears  separate  whenever  the  velocity  of  the  load  exceeds  that  of 

the  motor,  which  depends  upon  the  inertia  and  friction  distribution 
between  motor  and  load.  In  this  case  the  inei'tia  of  the  motor  is  equal 

to  the  inertia  of  the  load,  but  the  load  friction  is  much  greater  than 

the  motor  friction.  For  this  reason  the  gears  do  net  separate  until  the 

motor  has  sufficient  error  signal  in  the  opposite  direction  to  cause  the 

"*"A  study  is  being  carried  on  concurrently  to  include  the  conserva¬ 
tion  of  energy  and  the  elastic  impact  case  by  Lt.  N.  0.  Anderson  and 
Lt.  T.  W.  Luckett  at  the  U.  S.  Naval  Postgraduate  School. 
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Fiq.  3  Phase  Trajectory  for  Showing  Backlash  effects  (Example  l) 


motor  to  decelerate  faster  than  the  load.  As  the  error  signal  (which 
is  determined  by  the  load  position)  increases,  the  motor  velocity 
decreases  to  zero  and  begins  accelerating  in  the  opposite  direction. 

When  the  load  and  motor  positions  differ  by  an  amount  equal  to  the  back¬ 
lash  angle  as  shown  at  point  C  and  C ' ,  the  gears  make  contact  and  the 
lav/  of  momentum  is  applied  causing  the  motor  and  load  to  immediately 
assume  the  same  velocity,  that  velocity  being  dependent  upon  the  rela¬ 
tive  masses,  or  inertias,  of  the  motor  and  load.  This  new  velocity  is 

t 

shown  at  points  D  and  D  .  The  motor  and  load  now  both  move  in  accord¬ 
ance  with  the  combined  motion  equation  but  in  the  opposite  direction 
and  separated  by  an  amount  equal  to  the  magnitude  of  the  backlash,  as 
shovm  from  D*  to  S'  and  from  D  to  E. 

The  same  cycle  of  events  occurs  throughout  the  motion  of  the  system. 
The  effect  on  the  first  overshoot  is  observed  from  this  much  of  the  tra¬ 
jectory. 

Figure  4  shows  the  effect  of  backlash  on  a  system  which  has  a  load 
inertia  greater  than  the  motor  inertia  end  a  motor  friction  greater  than 
the  load  friction,  with  the  same  backlash  as  in  the  first  example.  It 
is  important  to  notice  at  point  B  that  the  gears  separate  earlier  than 
in  the  previous  example.  Both  the  heavier  load  inertia  and  lighter  load 
friction  contribute  to  its  slower  deceleration.  As  before,  when  the 
motor  and  load  are  separated  by  the  amount  of  the  backlash  the  gears 
make  contact  and  the  motor  and  load  assume  the  same  velocity  as  shown 
by  point  D  and  Dv.  Note  that  the  combined  velocity  is  positive  in 
this  case  because  of  the  high  load  momentum.  Thus  the  combined  motion 
is  in  the  positive  direction  with  the  error  signal  causing  the  motor  to 
act  as  a  brake  on  the  system  motion. 
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F i6. 4  Phase  trajectory  for  showing  backlash  effects  (example  Z) 


These  two  cases  serve  to  indicate  the  fact  that  every  case  nust  be 
studied  individually  since  every  combination  of  inertias  and  frictions 
cause  considerably  different  results. 

In  order  to  make  the  problem  mathematically  tractible,  the  follow¬ 
ing  assumptions  have  been  made: 

1.  Backlash  is  the  only  non-linearity  present  in  the  system. 

2.  Backlash  is  present  only  in  the  gears  coupling  the  motor  to 
the  load. 

3.  The  gear  contact  is  plastic  thus  applying  the  lav/  of  conserva¬ 
tion  of  momentum  only. 

It  is  believed  that  the  results  obtained  v/ith  these  assumptions 
give  a  good  approximation  of  the  effects  of  backlash  on  a  second  order 
servo  system. 
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Chapter  IV 

SOLUTION  OF  EQUATIONS  OF  MOTION 


When  operating  as  a  combined  system  the  servo  motion  can  be 


described  by  the  differential  equation. 


6L  +  2fo±  &i  +  *4*6L  =  oj/Ge 


(iv-l) 


When  the  load  drifts  free  within  the  backlash  region  the  load  motion 
follows  the  differential  equation. 


(IV-2) 


and  the  equation  for  the  motion  of  the  motor  is. 


/ 


(IV-3) 


The  characteristics  of  the  second  order  servo  with  backlash  is 
independent  of  the  undamped  natural  frequency.  The  proof  is  develop- 
ed  in  Ref.  2.  Therefore,  all  of  the  problems  solved  and  used  for  the 
construction  of  the  design  charts  are  based  on  an  undamped  natural  fre¬ 
quency  of  one  radian  per  second.  The  input  for  all  problems  was  a  step 
input  of  one  radian.  Thus  equation  (l)  becomes. 


(IV-11) 


+  oL  -  / 


The  differential  equations  IV-1^,  IV-2,  and  IV-3,  for  the  second 


order  servo  with  backlash  were  solved  on  the  Control  Data  Corporation 
Model  1604  Computer.  The  program  used  was  based  on  the  Runge-Kutta  Gill 
numerical  integration  method  of  solving  differential  equations.  A  com¬ 
plete  description  of  the  program  and  an  explanation  of  the  solution  are 
given  in  Ref.  2.  The  program  is  included  in  Appendix  C.  By  merely  insert- 
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ing  into  specific  memory  locations  the  system  parameters  for  the  problem 
to  be  solved,  the  solution  could  be  obtained  for  any  desired  second  order 
servo. 

A  study  of  the  describing  differential  equations  reveals  that  there 
are  five  system  parameters  which  are  independent.  These  are  F^,  J. , 
F^»  and  A  .  To  facilitate  the  graphing  of  the  results,  ratios  of 
J  F 

L  and  L  were  selected,  and  from  these  selected  ratios  discrete 
JM  Ftotal 

values  of  the  parameters  J^,  F^,  J^, ,  and  F^,  were  calculated  for  use  in 
the  computer  solution. 

The  system  damping  coefficient  is  determined  by  the  relationship, 

'~'T 

Taking  J,p  =  =  1  simplified  calculations,  but  did  not 

destroy  the  generality.  The  system  problem  was  solved  for  all  combina¬ 
tions  of  the  following  ranges  of  inertia  and  friction  ratios  and  damping 
coefficient; 


£l  =  1/50,  1/9,  1/4,  1,  4,  9,  and  90 

JM 

M 

fl 

~~  s  0.1,  0.2,  0.4,  0.6,  0.8,  1.0 

* total 

f  =  0.1,  0.2,  0.3,  0.4,  0.5,  0.6,  0.8,  1.0 


The  manner  in  which  the  results  are  presented  in  the  final  design 
charts  permits  interpolation  through  the  practical  range  of  system 
parameters. 

The  amount  of  backlash,  A  ,  used  in  the  majority  of  the  problems 
was  0.3  radian  and  it  was  a  purpose  of  this  thesis  to  determine  the 
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effect  of  other  magnitudes  of  backlash.  At  the  outset  several  problems 
were  solved  varying  the  amount  of  backlash  while  holding  the  other  para¬ 
meters  at  constant  values.  Figures  5-7  are  plots  of  the  results  of  these 
solutions.  Two  items  of  significant  interest  were  noted  from  these  re¬ 
sults,  First,  the  amplitude  of  maximum  overshoot  varies  linearly  with 
the  magnitude  of  the  backlash  angle.  This  fact  is  significant  since  it 
permitted  the  use  of  a  constant  backlash  angle  in  the  solution  of  the 
problems,  yet,  in  effect,  produced  results  for  any  amount  of  backlash. 
Secondly,  the  amplitude  of  maximum  overshoot  at  zero  backlash  angle 
agrees  with  that  obtained  by  the  analytical  solution  of  the  second 
order  linear  system.  This  agreement  is  taken  as  a  verification  of  the 
accuracy  of  the  computer  solutions. 
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Chapter  V 


DISCUSSION  OF  RESPONSE  CHARACTERISTICS 

A.  Maximum  Overshoot  for  the  Plastic  Impact  Case. 

After  numerous  data  had  been  accumulated,  several  schemes  of  pre¬ 
senting  the  information  were  attempted  to  find  an  optimum  method  for 
use.  One  such  attempt  is  shown  in  figure  8  which  is  a  plot  of  maximum 
overshoot  versus  the  ratio  of  load  friction  to  the  total  friction  at  a 
constant  damping  coefficient.  This  method  was  found  to  have  serious 
disadvantages.  Because  of  appreciable  crossing  of  the  lines  of  con¬ 
stant  inertia  ratio,  it  was  found  difficult  to  interpolate  throughout 
the  desired  range  of  inertia  ratios.  Also  there  was  generally  a  small 
variation  of  maximum  overshoot  with  friction  ratios. 

Another  method  which  was  attempted  was  a  plot  of  maximum  over¬ 
shoot  versus  damping  coefficient  at  constant  friction  and  inertia 
ratios.  These  curves  showed  the  effect  of  various  inertia  and  fric¬ 
tion  ratios  in  a  system  with  backlash  on  the  linear  system  curve  of 
overshoot  versus  damping  coefficient  (figure  9),  but  was  found  diffi¬ 
cult  to  use.  The  number  of  curves  needed  to  present  the  data  would 
have  been  excessive.  Also  interpolation  throughout  the  desired  range 
of  inertia  ratios  would  have  been  extremely  difficult. 

The  curves  of  figures  10  through  15,  which  are  plots  of  maximum 

overshoot  versus  were  found  to  be  the  most  desirable  method  of 

L  M' 

presenting  the  data.  None  of  the  disadvantages  apparent  in  the  other 
methods  is  present  here.  From  these  six  charts  an  accurate  prediction 
of  the  maximum  overshoot  can  be  made  for  any  arbitrary  second  order 
servo  system  with  backlash  using  the  following  equation; 


19 


M 


pt 


M 


|A=A, 

'f=// 


pt 


A-<9 


A/o.3 


H 


M 


(V-l) 


ptU*3  ptU =0 

The  use  of  these  charts  in  predicting  system  response  is  explained  and 

'V 

\{ 

and  thi 


illustrated  in  Appendix  A, 

As  was  pointed  out  earlier,  the  solutions  were  made  and  th£'  results 

were  plotted  using  a  backlash  angle  of  0.3  radian.  Thus  the  values  of 

amplitude  of  maximum  overshoot  on  the  curves  are  significant  only  in  that 

they  represent  the  slope  of  the  overshoot  versus  backlash  straight  lines. 

This  slope  must  be  calculated  in  order  to  scale  the  value  of  overshoot 

read  from  the  curve  for  the  magnitude  of  backlash  actually  expected. 

Therefore,  in  order  to  refine  the  prediction  scheme,  other  design  curves, 

figures  16  through  21,  were  calculated  using  the  curves  of  figures  10 

clM 

through  15  just  discussed.  Along  the  curves  the  slope  pt 

d  A 

calculated  by. 
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These  values  of  dM  .  /da.  xvere  then  plotted  versus  inertia  ratios  at 

pt' 

constant  friction  ratios  as  before.  These  curves  give  a  more  straight¬ 
forward  prediction  method  with  the  following  equations 
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The  use  of  these  design  charts  is  explained  and  illustrated  in  Appendix  A 
along  with  the  former  method. 

Figure  23  shows  a  graphical  presentation  of  the  transient  response 
of  the  system  containing  various  parameters.  A  comparison  of  curves  (3) 
and  (4)  readily  shows  the  effect  of  backlash  on  the  maximum  overshoot 
with  ether  parameters  remaining  constant.  A  comparison  of  curves  (1) 
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and  (4)  shews  the  effect  of  varying  the  damping  ratio  of  the  system  with 
other  pax  ameters  remaining  constant. 
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B.  Maximum  Overshoot  for  the  Elastic  Impact  Case, 

One  of  the  assumptions  made  in  the  beginning  of  this  study  was 
that  the  contact  of  the  gear  teeth  is  plastic  with  no  bouncing  of  the 
gears o  Since  this  study  was  begun,  an  investigation  has  been  completed 
at  the  Uo  3.  Postgraduate  School  which  considers  the  backlash  case  with 
varying  amounts  of  resilience  in  the  gear  teeth.  This  investigation, 
carried  on  by  T„  W.  Luckett ,  and  N.  0.  Anderson  has  produced  some  results 
bearing  on  maximum  overshoot  which  will  be  included  here.  The  informa¬ 
tion  ,  while  not  included  as  part  of  their  study,  was  available  from  the 
data  taken  from  the  computer  solution  of  their  problem. 

The  analysis  of  this  data  revealed  a  surprising  effect  of  gear 
bounce  on  the  maximum  overshoot  of  the  system.  The  amount  of  over¬ 
shoot  with  plastic  impact  is  always  greater  than  that  with  elastic 
impact.  The  difference  between  this  maximum  overshoot  decreases  as 

F,/F,  ,  .  increases.  However,  under  no  circumstances  was  the  effect 

U  total  5 

on  maximum  overshoot  greater  than  .01  radians.  The  effect  on  the 
time  of  max  mum  overshoot  appears  negligible.  In  a  few  cases  the  time 
was  different  for  the  plastic  and  the  elastic  cases  but  the  difference 
was  very  small  and  no  pattern  of  difference  could  be  ascertained. 

These  re  stilts  point  out  the  validity  of  the  original  assumption 
made  for  this  study.  In  a  practical  system  there  will  exist,  of 
course,  some  amount  of  gear  bounce  when  the  gears  contact  each  other. 

The  effect  cr»  maximum  overshoot ,  however,  is  so  small  that  the  curves 
obtained  in  this  study  remain  reasonably  accurate. 
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G o  Settling  Time  ana  Time  of  Maximum  Overshoot . 

1,  Settling  Time 

In  determining  the  effect  of  backlash  on  the  settling  time  of  a 
second  order  servo  system,  it  was  necessary  to  first  find  a  definition 
for  settling  time,  which  would  be  suitable  for  all  of  the  cases  which 
•were  investigated.  Normally,  settling  time  is  defined  as  the  time  for 
the  output  to  reach  98%  of  its  final  value.  This  definition  has  real 
meaning,  however,  only  for  linear  stable  systems  in  which  there  are  no 
sustained  oscillations.  For  systems  which  do  exibit  small  sustained 
oscillations,  the  settling  time  is  sometimes  defined  as  the  time 
required  for  the  amplitude  of  oscillations  to  decrease  to  within  2% 
of  the  final  value.  For  very  small  amplitudes  of  oscillations  this 
definition  gives  an  unrealistic  value  for  settling  time.  For  the 
purpose  of  this  study  settling  time  is  considered  as  the  time  required 
for  the  system  output  to  settle  to  98%  of  the  difference  between  the 
magnitude  of  the  input  and  the  amplitude  of  the  sustained  oscillations 
as  shown  in  the  sketch  below. 


3? 


For  all  cases  for  which  the  damping  coefficient  was  less  than  0,6 
there  were  no  conclusions  as  to  effects  of  backlash  on  the  system 
settling  time  since  the  magnitudes  of  the  sustained  oscillations  were 
always  quite  large  for  the  amount  of  backlash  used  in  the  solutions  and 
the  settling  time  was  difficult  to  obtain  accurately.  For  those  cases 
for  which  the  damping  coefficient  was  greater  than  0.6*  the  sustained 
oscillations  were  small  if  they  existed  and  an  analysis  revealed  the 
following  general  conclusions s 

1.  As  the  size  of  the  backlash  angle  increases*  the  settling  time 
will  increase  provided  the  other  parameters  are  kept  constant.  For  the 
cases  investigated*  this  increase  in  settling  time  was  nearly  proportion¬ 
al  to  the  backlash  angle. 

2.  As  the  ratio  J^/J^  increases*  with  other  parameters  constant* 
the  settling  time  increases. 

3.  At  the  ratio  F, 

constant*  the  settling  time  decreases. 

The  data  included  in  Tables  I-III  show  these  general  conclusions 
and  also  some  quantitative  effects  of  backlash  on  settling  time. 


/Ftotalis 


increased*  with  other  parameters 


38 


TABLE  I 


SETTLING  TIME  DATA 


f 

Jr 

Vj 

M 

fl/ 

Ftotal 

Radians 

A 

Seconds 

t 

s 

.6 

1/9 

.2 

.3 

28 

.6 

1/9 

•  4 

.3 

21 

.6 

1/9 

.6 

.3 

16 

.6 

1/9 

.8 

.3 

13 

.6 

1/9 

1.0 

.3 

7 

o  w 

1/4 

.2 

.3 

28 

.6 

1/4 

.4 

.3 

24 

.6 

1/4 

o 

.3 

16 

.6 

1/4 

.8 

»3 

14 

.6 

1/4 

l.C 

.3 

9 

.6 

1 

„2 

.3 

29 

.6 

1 

.4 

.3 

28 

.6 

1 

.6 

.3 

23 

.6 

1 

.8 

.3 

17 

.6 

1 

l.C 

.3 

10 

.6 

4 

.2 

.3 

29 

.6 

4 

.4 

.3 

28 

.6 

4 

.6 

.3 

26 

.6 

4 

.8 

»3 

23 

.6 

4 

1.0 

.3 

10 

.6 

9 

,2 

.3 

29 

.6 

9 

•  4 

.3 

29 

.6 

9 

.6 

.3 

26 

.6 

9 

.8 

.3 

23 

.6 

9 

1.0 

.3 

10 
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TABLE  II 


SETTLING  TIME  DATA 


f 

JL 

M 

F 

L, 

/  F 

total 

Radians 

A 

Seconds 

t 

s 

.6 

1/9 

.8 

.03 

5 

.6 

1/9 

.8 

o 

o 

7 

06 

1/9 

.8 

.15 

8 

.6 

1/9 

.8 

.30 

13 

.6 

1 

.8 

.03 

11 

.6 

.8 

.10 

14 

.6 

1. 

.8 

.15 

16 

.6 

1 

.8 

o 

cn 

o 

17 

.6 

4 

.8 

.03 

12 

.6 

4 

.8 

.10 

14 

.6 

4 

.8 

.15 

17 

o6 

4 

.8 

.30 

23 

.6 

9 

.8 

.03 

13 

.6 

9 

.8 

.10 

16 

.6 

9 

.8 

.15 

18 

.6 

9 

p 

fr  'w' 

.30 

23 
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TABLE  III 


SETTLING  TIME  DATA 


T 

% 

/jM 

fl 

L/ 

•  F 

total 

Radians 

A 

Seconds 

.8 

1/9 

.2 

.3 

8 

.8 

1/9 

.4 

*3 

7 

.3 

1/9 

.6 

.3 

6 

.3 

1/9 

.8 

.3 

5 

.8 

1/9 

1.0 

.3 

5 

.8 

1/4 

.2 

»3 

20 

.8 

1/4 

.4 

.3 

9 

08 

1/4 

.6 

.3 

7 

08 

1/4 

.8 

.3 

6 

.8 

1/4 

1.0 

.3 

5 

.8 

1 

.2 

.3 

43 

.8 

1 

•  4 

.3 

31 

.8 

1 

.6 

.3 

10 

.8 

1 

.8 

.3 

8 

o  8 

1 

1.0 

.3 

6 

.8 

4 

.2 

.3 

48 

.8 

4 

.4 

.3 

36 

.8 

4 

.6 

.3 

12 

.8 

4 

.8 

.3 

9 

.8 

4 

1.0 

.3 

6 

c  8 

9 

.2 

.3 

52 

.8 

9 

.4 

.3 

43 

.8 

9 

.6 

.3 

13 

08 

9 

.8 

.3 

9 

.8 

9 

1.0 

.3 

6 
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2.  Time  of  Maximum  Overshoot 


For  each  of  the  cases  investigated  the  time  of  maximum  overshoot 
(peak  time)  was  recorded  from  the  printout  of  the  solution.  It  was 
found  that  the  presence  of  backlash  causes  a  slight  increase  in  the 
peak  time  over  the  linear  system  for  those  cases  where  and 

has  negligible  effect  in  those  where  'S  K.  .5  .  When  F^/F+  is 
increased  from  zero  to  0.5*  other  parameters  constant  and  J*  .5* 
peak  time  increases  slightly  and  remains  essentially  constant  for 
Fl/Ftotal  ^  ‘  As  Jl/JM  is  increased,  other  parameters  constant 

and  'S  y  .5*  peak  time  is  decreased. 
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D.  Other  System  Characteristics . 

A  careful  study  of  the  results  as  plotted  in  figures  10-15  reveals 
some  interesting  phenomena.  At  all  friction  ratios  it  is  observed  that 
as  the  ratio  of  load  inertia  to  motor  inertia  is  made  smaller  the  ampli¬ 
tude  of  the  overshoot  is  less.  There  is  a  point  where  the  overshoot  is 
identical  to  that  of  a  linear  second  order  system*  in  other  vrords,  where 
backlash  has  no  effect  on  the  overshoot.  Note  that  this  point  occurs 
at  different  inertia  ratios  for  each  value  of  damping  ratio.  Also  this 
point  varies  depending  on  the  value  of  friction  ratio.  If  the  inertia 
ratio  is  decreased  beyond  this  point  then  the  amplitude  of  overshoot 
remains  constant  at  the  value  of  the  linear  case.  Figure  22  is  a  plot 
of  the  boundaries  of  the  combinations  of  variables  which  will  cause  an 
effect  on  the  amplitude  of  overshoot  if  the  system  has  any  appreciable 
amount  of  backlash.  The  coordinates  are  friction  ratio  and  inertia 
ratio  with  the  boundaries  labeled  with  values  of  damping  ratio.  Areas 
to  the  right  of  a  boundary  are  combinations  of  the  variables  which 
affect  overshoot.  These  curves  were  plotted  using  information  deducted 
from  figures  10-15.  Three  points  were  selected  as  test  points  to  prove 
the  validity  of  the  boundaries  as  plotted.  These  points  are  shown  on 
the  figure  as  points  "A",  "B",  and  "C".  The  system  parameters  used  at 
the  test  points  were; 


Point  A  Point  B  Point  C 


f 

=  0.S 

f 

=  0.6 

f 

-  0.3 

Fl/Ftotal 

=  O.S 

Fl/Ftotal 

=  0.6 

FL/Ftotal 

=  0.4 

JL/JM 

=  1.0 

VJM 

-  0,25 

VJM 

=  0.111 
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The  system  was  solved  on  the  digital  computer  to  obtain  the  maximum 
overshoot  for  these  three  sets  of  parameters.  According  to  the  curves 
point  A  falls  to  the  left  of  the  boundary  of  the  f  =  ,S  area. 
Therefore,  the  overshoot  should  be  the  same  as  that  of  the  linear  case. 
The  computer  solution  gave  an  overshoot  of  1.019,  which  is,  in  fact, 
identical  to  that  of  the  linear  case  as  seen  in  figure  9.  Similarly, 
point  B  falls  to  the  right  of  the  f  -  0.6  boundary  and  in  the  area 
where  the  amplitude  of  overshoot  should  not  be  the  same  as  for  the 
linear  case.  The  computer  solution  gave  a  maximum  overshoot  of  1.110 
which  is  slightly  greater  than  the  linear  case.  Point  0  falls  to  the 
right  of  the  f=  0.3  boundary  and  in  the  area  where  the  amplitude  of 
overshoot  should  be  different  from  that  of  the  linear  system.  The 
computer  solution  gave  a  maximum  overshoot  of  1.41  while  that  of  the 
linear  system  is  1.372.  From  these  boundary  curves  one  cannot  deter¬ 
mine  the  exact  effect  of  backlash  on  the  maximum  overshoot  of  a  system, 
but  can  determine  in  which  systems  backlash  will  cause  an  appreciable 
effect. 

Another  phenomena  which  is  revealed  in  the  results  as  plotted  in 
figures  10-15  is  that  for  all  cases  the  greatest  maximum  overshoot 
occurs  in  the  range  0.1  J^/J^  ^  10. 0  for  the  range  of  J^/J^  test¬ 

ed.  However,  it  appears  that  since  the  curves  are  trending  upwards  at 
the  high  end  of  the  axis  that  at  very  large  inertia  ratios  the 

overshoot  would  be  greater.  It  is  felt  that  these  large  inertia  ratios 
are  outside  the  practical  range.  An  attempt  was  made  to  explain  these 
phenomena  but  without  success.  It  is  often  very  difficult,  if  not 
impossible, to  find  a  physical  explanation  for  certain  behavior  caused 
by  non-linearities.  They  can  merely  be  pointed  out  from  analyses  of 


experimental  information  and  accounted  for  accordingly. 

It  is  interesting  to  note  that  when  the  ratio  of  load  friction  to 
motor  friction  is  below  0.4,  as  in  figures  10  and  11,  there  is  a  possi¬ 
bility  of  an  overshoot  for  the  critically  damped  case,  i.e.  "S  =  1.0, 
depending  on  the  inertia  ratio.  At  first  glance  one  would  reason  that 
if  the  system  is  critically  damped  there  should  be  no  overshoot.  How¬ 
ever,  it  can  be  seen  that  if  the  friction  ratio  is  small  and  the  inertia 
ratio  is  large  then  the  ratio  of  which  is  the  slope  of  the  drift¬ 

ing  load  trajectory,  is  made  small.  Therefore,  when  the  load  separates 
from  the  motor,  its  velocity  remains  fairly  constant  at  a  relatively 
high  value  and  its  momentum  carries  it  over  the  input  to  cause  an 
overshoot. 
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Chapter  VI 
CONCLUSIONS 

The  following  is  a  summary  of  facts  which  can  be  concluded  from  the 
results  of  this  investigation: 

1.  The  amplitude  of  maximum  overshoot  in  a  second  order  servo  sys¬ 
tem  with  backlash  bears  a  linear  relationship  with  the  magnitude  of  the 
backlash  angle.  The  slope  of  the  straight  line  relating  overshoot  to 
the  backlash  angle  varies  depending  on  the  division  of  inertia  between 
motor  and  load,  the  division  of  friction  between  motor  and  load,  and  the 
system  damping  coefficient.  The  design  charts  included  in  this  thesis, 
drawn  from  the  results  of  many  solutions  of  the  equations  of  the  motion 
of  the  system  and  based  on  a  backlash  angle  of  0.3  radian,  can  be  used 
to  predict  the  amount  of  overshoot  of  a  particular  servo  system. 

2.  The  fact  that  the  gear  teeth  have  a  certain  amount  of  resilience 
and  the  contact  is  not  plastic,  has  negligible  effect  on  the  results  ob¬ 
tained  in  this  study. 

3.  The  primary  effect  of  backlash  on  the  settling  time  of  a  second 
order  servo  is  to  increase  it,  the  increase  being  proportional  to  the 
magnitude  of  the  backlash  angle.  Due  to  the  presence  of  limit  cycles 

of  sizable  magnitudes  for  systems  with  -J  < .6,  settling  time  for  those 
systems  was  difficult  to  define  and 'no  conclusions  were  reached.  For 
those  systems  with  (a)  Increasing  F^/F^^^  decreases  settling 

time,  (b)  Increasing  J./J^  increases  settling  time,  (c)  Increasing  the 
damping  coefficient  increases  settling  time,  and  (d)  Increasing  back¬ 
lash  increases  settling  time. 

4.  The  presence  of  backlash  in  a  system  causes  a  slight  increase 
in  the  time  of  maximum  overshoot  (peak  time)  for  those  systems  whereby 
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and  has  negligible  effect  in  those  where  When  Ft/F.  ,  is 

L  total 

increased  from  zero  to  0.5,  other  parameters  constant  and  jr>.5, 
peak  time  increases  slightly  and  remains  essentially  constant  for 
F^/F^+  ^  >  .5.  As  is  increased,  other  parameters  constant  and 

f  >  .5,  peak  time  is  decreased. 

5.  With  certain  combinations  of  parameters  present  in  a  system, 
backlash  effects  no  change  in  the  amplitude  of  overshoot  over  that  of 
the  linear  system.  Regions  of  these  combinations  are  presented  in  chart 
form  in  figure  22. 

6.  For  all  friction  ratios  and  damping  coefficients  the  greatest 
maximum  overshoot,  when  backlash  is  present,  occurs  in  the  range 

0.1  <  10.0  for  the  range  of  J^/J^  investigated. 

7.  In  systems  where  ^jA’total  ^  w*ien  ba°klash  is  present, 

there  will  be  an  overshoot  even  for  critically  damped  systems  if  J^/J^ 


is  high  enough. 
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APPENDIX  A 


USE  0?  DESIGN  CHARTS 


In  order  to  demonstrate  the  procedure  for  using  the  design  charts 
shovm  in  figures  10-21,  an  example  problem  is  solved  obtaining  the 
maximum  overshoot,  M  for  a  second  order  servo  with  backlash.  From 
Ref.  2,  the  differential  equation  for  the  combined  motion  of  an  armature 
controlled  motor  is: 

e  *  <*■« 


Given  the  following  specifications  of  a  D.  C.  shunt  motor  coupled 


through  a  gear  box  to  a  load: 


Horsepower 
Rated  RFM 


R 

K 

i 

K. 


m 


m 


K 


/° 

A 


=  .008 
=  4000 
=  282  ohms 

=  .233  volts/radian/second 

=  .1464  lb. -feet/ampere 
=  8.2  X  10“6  slug-feet2 
=  3.34  X  10“^  lb.-feet/radian/sec. 
=  10  volts/radian 
=  1 
=  .05 

=  16.4  X  10“^  slug-feet2 
=  10.4  X  10~6  lb. -feet/radian/sec. 


52 


The  required  values  shown  in  Equation  A-l  are  obtained  as  follows? 


Ke  =  ■  (±64  *  / 'O 

/°  ^  0)x2  #2 

iZ,  _  f.2  Jr  /O'6 

/>*  OJ* 


0052  18- FT 

f-3  SL  L/<S  -  FT 


COn 


.OOS2 


f-2  *  /O'6  y  /6. 4  */0~& 


/  4.5 


/Q£\P 

Sec. 
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7° 


/o.  4  x  /O  y- 
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(O' 


7totai  “  ft  +  Fi  ~  {^3.74  ^  /2/)  x  /o'6  =  /S5  *  /o'6  5±£L 

'CM/se-c 

-6  S.3-XT 


kt  45  _  .  (464  4  >j233. 

0)2  x  2  S2 


f  -  /FL-EL .  /_L  )  =  x_9_ 
+  l_Z  ^ 


=  A?/  /♦* 


^4  D/s £c 


/6.4  2  /O 


JL/JM  /2  x  /o~ ' 


-  6 


=  o 


fl/f 


tot  3.1 


'3S  x  /o  ~6 


•  <?77 


Two  methods  are  shown  for  obtaining  the  maximum  overshoot  with  a 
backlash  angle  of  .05  radians  in  the  gear  trains. 

The  first  method  is  accomplished  using  the  curves  of  figures  10-15. 
The  procedure  is  simple,  but  requires  more  steps  than  does  the  second. 
Also  one  must  remember  that  these  curves  are  drawn  for  a  backlash  of 
0.3  radian  and  must  be  scaled  for  the  expected  amount  of  backlash.  The 
second  method  makes  use  of  the  design  charts  of  figures  16-21, 
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Method  1 


Entering  figure  10,  (^/^total  =  0,1^#  at  di/dM  ~  2o0>  and  a  damp‘ 
ing  ratio, S  =  0.19,  the  amount  of  overshoot  for  a  backlash  angle,  -A  - 
0.3  rad.  is  M  ,  =  1.653. 

pt 

Using  Equation  V-l 


=  M  , 
pt 

A=.05 

+  .05/.3 

M  .  -  M  , 

A-0 

f^.!9 

"f  *  .19 

\fsJ9 

f=J$ 

From  the  second  order  linear  characteristics  curve  of  figure  9 


Mpt  =  1.545 

hence, 

M  .  =  1.545  +  .05/.3  1.658  -  1.545 

pt. 

=  1.545  +  .0189 

-  h&k 

Method  2 

Using  Equation  V-3 


1  A=.05 

M  . 
pt 

+  dM  . /d A 
pt 

\fz./9 

X 

f=./$ 


A 


from  figure  16,  (FL/Ftotal  =  °*1^  at  Ji/JM  =  2*°>  ^  a  damp±ng  ratio> 


=  0.19,  we  obtain  dM^/dA  =  0.357 


From  figure  9 

M  .  =  1.545 

pt  I 

I  f=  .  19 

hence, 

Mpt  -  1.545  +  .357  X  .05 

=  1.545  +  .018 

=  1.563 

This  compares  closely  to  the  answer  obtained  using  Method  1. 
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The  procedure  for  using  the  charts  remains  unchanged  for  servo 
systems  which  utilize  other  types  of  motors.  The  only  difference  which 
arises  is  in  the  calculation  of  the  chart  parameters  from  the  known 
system  constants.  For  a  field  control  motor  if  the  field  time  constant 
is  negligible,  the  equation  is  similar  to  the  armature  control  motor 
equation. 


9  f 


E 

Em  -f- E 


G  f 


— 9  = 


(A-2) 


In  like  manner  the  equation  for  the  system  with  a  two  phase  motor  is 


Q 


-h 


ElL* 


6 


+  — 9 


L  T  (A-3) 

Em  Ei. 


where 


K 


n 


K 


m 


dT 
~  bn 

EL 

2>e, 


(from  characteristic  curves  of  two  phase  motor) 
(from  characteristic  curves  of  two  phase  motor) 


The  calculation  of  the  chart  parameters  is  a  simple  matter  once  the 
defining  equations  are  known . 
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